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Abstract 

The BFKL equation and the /cT-factorization theorem are used to obtain predictions for F2 
in the small Bj0rken-x region over a wide range of Q^. The dependence on the parameters, 
especially on those concerning the infrared region, is discussed. After a background fit 
to recent experimental data obtained at HERA and at Fermilab (E665 experiment), we 
find that the predicted, almost independent BFKL slope A > 0.5 appears to be too 
steep at lower values. Thus there seems to be a chance that future HERA data can 
distinguish between pure BFKL and conventional field theoretic renormalization group 
approaches. 



1 Introduction 



The usual kinematic variables used for discussing deep-inelastic electron-proton scattering 
are derived from the four-momenta p of the incoming proton and q of the exchanged virtual 
photon: = —q^ and the Bj0rken variable x = Q"^ /2p-q. In the region where higher twist 
effects are likely to be negligible, i. e., for W'^ = Q"^ 1) ^ 10 GeV^, the (Dokshitzer- 

Gribov-Lipatov-)Altarelli-Parisi [(DGL)AP] set of equations [0, describes the evolution 
of the structure function F2 with very well. In leading order (LO) all powers of 
Us ln((5^//i^) are summed by the AP evolution equations, which take into account just 
strongly ordered parton-fej- ladders. Nowadays usually the next-to-leading order (NLO) 
set is used, where terms of the form a" ln"~^((5^/ /i^) are also summed, taking non-ordered 
contributions (covariantly) into account as well. Even in NLO, the AP equations do 
obviously not contain all leading logarithms in x. Thus one might naively expect the AP 
framework to break down at some small value of x, where a resummation of all powers 
of as ln(l/x) should be necessary, although no perturbative instability between LO and 
NLO has been observed thus far in the presently relevant kinematic regime 0, ^, 0. 

Such a resummation in LO is provided by the Balitskij-Fadin-Kuraev-Lipatov (BFKL) 
equation . The equation treats the gluons only, which are expected to be the dominant 
partons at small x. This might be deduced from the AP splitting functions which 
become ~ 1/x for splitting to gluons and constant for splitting to quarks in the small 
X limit. It should be remarked though, that this argument may be too simple, as it 
neglects the influence of the particular shape of the parton distributions used |^ and 
furthermore does not respect the fundamental energy- momentum conservation constraint. 
The BFKL resummation is formally correct for all Q^, but a fixed coupling was used in 
its derivation. Furthermore it is based on LO perturbative QCD, using (non-covariant) 
cut-off regularizations, thus we are limited to sufficiently large . Since the resummation 
assumes that subleading terms in x are small, including those involving logarithms of Q^, 
it will also not be valid at high . A conservative range for would be 4 to 50 GeV^, 
and the range 0.8 to 120 GeV^ explored in this work should be regarded as the extreme 
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limit. 



To estimate the effects of tlie subleading terms, a NLO resummation in x would be 
necessary. But even though there has been some progress in that direction 0, a final 
result has not yet been obtained. Ultimatively it should even be possible to find a unified 
evolution equation covering the whole perturbative region [|, [1^]. But a calculation 
that can be confronted with experiment is still missing, thus we stay with the usual BFKL 
formalism for the time being. 

Since a fixed coupling constant seems unreasonable in view of the running coupling 
of the AP framework we wish to connect to, the replacement tts(Q^) is done by 

hand. There is really no rigorous motivation for this step. Some trust in this procedure 
can be gained by considering the representation of the evolution equations as ladder 
diagrams. It is well known that the LO AP equations can be represented in a physical 
gauge by a sum of ladder graphs which are strongly ordered in the transverse momentum 
k'^: > /c^ > > . . . > /^i > ^0- small x limit of the AP equation 

we consider only the dominant gluon ladder (see Fig. |l]) and keep only the terms with 
double leading logarithms (DLL) as ln(l/a;) ln(Q^/yU^). This corresponds to introducing 
an additional strong ordering in x: x ^ Xn-i ^ . . . ^ Xq. The BFKL evolution can also 
be described in terms of a (reggeized) gluon ladder, using the strong ordering in x only, 
to get the leading logarithms in x. If we let the coupling run, then BFKL will reduce to 
LO DLL upon imposing the ordering in transverse momentum. 

The main feature of the BFKL evolution with fixed coupling constant is the growth 
of the unintegrated gluon ~ x~^ with A = ln2 ^ 0.5. Due to the dominance of the 
gluons at small x we expect a corresponding rise at small x in the structure functions 
too, once the off-shell gluons have been appropriately coupled to the quark sector. This 
expectation has been confirmed [|ll|, [1^ for a small range of medium even in the case 
of running coupling BFKL, with F2 ~ x^^ and A > 0.5 for x < 10~^. 

Experimentally the situation has improved drastically since the advent of the HERA ep 
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collider. Before HERA only the Fermilab experiment E665 |T3| was able to reach the small 
X region, but at rather low Q^. Now HERA takes data [|14|, |l^ with x/Q^ >3 ■ 10~^ over 
a wide range of Q^. The observed strong rise of F2 at small x has boosted the interest in 
the BFKL formalism. On the other hand the dynamically generated AP partons |Q create 
a steep gluon differently, via a long evolution length in Q^, and successful parameter-free 
predictions [0, ^ |l^ have been given long before HERA started to operate. Alternatively, 
the present data can be fitted using the NLO AP evolution equations: Then a term 
~ x'^ has to be assumed for the gluon distribution, e. g. in MRSG [0. Possibly this 
term mimicks the BFKL behaviour. But this is not clear, since both methods describe 
the data equally well. 

We conclude that a detailed comparison of the standard BFKL formalism with the 
new data is necessary. Only this can tell us if BFKL can rival conventional field theoretic 
renormalization group (AP) evolution equations in describing the measured structure 
function F2. Our calculations are based on the methods employed by Askew, Kwiecihski, 
Martin and Sutton (AKMS) |]ll], |T^, which will be described briefly in the following. 



1.1 BFKL equation and /c^-factorization 

The unintegrated gluon distribution /(x, A;^), which is related to the familiar integrated 
gluon distribution used in the AP equations by 



n dp 

xg{x,Q^) = J —f{x,k^), (la) 



f{x,e) 



dxg{x, Q"^) 



(lb) 



ain(Q2 

depends on the transverse momentum k. Using it, one can write the BFKL equation as 

g 

/(x,fc2) = /o(x,P)+ f I dk"K{e,k'')f{y,k'), 

J X y J 



One could use a suitable input /o, the so called "driving term", and solve the equation 
iteratively [|10|, but this procedure allows no simple connection to the known AP region. 



Instead we can obtain an evolution equation in x from the integral equation by differ- 
entiating with respect to ln(l/x) 



X = k 



dx Tx I k'"^ 



\k'^-k^\ {Ak'^ + k^)-2 



(2) 



assuming that the derivative of the driving term /o can be neglected. Since the driving 
term describes the gluon content without any BFKL evolution, it is reasonable to assume 
that it is connected to the non-perturbative "soft" pomeron [^. Because of the soft 
pomeron's weak x dependence ~ we expect df^/dlnil/x) to be small. 

Eq. (H) can then be used to evolve the unintegrated gluon to smaller x, using a suitably 
modified AP input as boundary condition at Xq = 10^^ by applying Eq. ([Tap. An obvious 
problem in Eq. (|^) is posed by the integration over k'"^, which starts at zero. Even for 
the AP gluon distributions we use in this work the limit of validity is about 1 GeV^. We 
employ a simple ansatz to continue the gluon into the infrared (IR) region, which will 
be described later in this paper. A similar comment applies to the upper limit of the fc'^ 



integration. The upper limit introduced by energy conservation is close to infinity |Tl|, |T^ , 
so that for practical calculations an artificial ultraviolet cutoff has to be introduced. The 
dependence on the IR and UV treatment will be thoroughly discussed later on. 

To obtain predictions for F2, we have to convolute the (off-shell) BFKL gluon, i. e. 
the gluon ladder in Fig. with the photon-gluon fusion quark box F^^^ using the kx- 
factorization theorem P, 



dk'^ dy 



X 



F^M) = J ^ 7^? ; ^^^°^(2/,A:^g^), (3) 

with i = T,L denoting the transverse and longitudinal parts, respectively. The expressions 
for F^^^ can be found in and references therein. The BFKL prediction for F2 is then 



simply the sum of the calculated Ft and F^. Obviously here the same problems with the 
/c'^ integration occur, which are circumvented by the methods mentioned above. 

It is important to notice that the BFKL gluons are off-shell (A;^ 7^ 0). The term F^'^'^ / k"^ 
in the above equation then corresponds to the structure function of a virtual gluon. In 
contrast the AP formalism is based on on-shell gluons, which is a good approximation 
due to the strong ordering in transverse momentum encountered in AP evolutions in LO. 
This strong ordering allows us to perform the k'^ integration in Eq. (j^). Thus, ignoring 
complications due to the coUinear singularities for simplicity, one arrives at the usual mass 
factorization equation 

F^{x,Q')= f ^g(-,QA F,{y,Q'). 

Jx y \y J 

Here Fi plays the role of the on-shell gluon structure function, whose x dependence stems 
from the AP splitting function Pqg, and g is the integrated gluon. In NLO AP evolutions 
the strong ordering in k"^ does not hold anymore due to the emission of a second gluon, 
but the interacting gluon is still considered to be on-shell in comparison with the hard 
scattering scale Q^. 



This shows that it is inconsistent to simply feed the evolved BFKL gluon via Eq. ([Ta|) 
back into the AP equations below the limit set by Xq. Calculations attempting to use 



BFKL gluons below and AP gluons above xq within the AP formalism ignore the 



essential ojj-shellness of the BFKL gluons. This casts first doubts on a recent BFKL 



analysis of the HERA HI data using this method |^ . No such problem persists when we 
use just this gluon mixture to drive the general fcr-factorization Eq. (|^), which reduces 
to the mass factorization in the AP region. 

Even though the dominant contribution at small x should come from the BFKL gluons, 
a certain amount of "background" in F2 due to quarks and non-perturbative effects should 
be taken into account. We expect the background to be comparatively small and also to 
vary much less with x. We shall use an ansatz motivated by the soft pomeron C/px~°'°^, 
where the constant Cjp is fitted to the data. After this general outline of our method, we 



will now proceed to a detailed discussion of the underlying formalism. 



2 Suitable input for the BFKL evolution 



We focus our analysis on the gluon distribution used in ||12|, i.e. a gluon based on the 



MRS Z^o-set of parton distributions |2^, but evolved with the leading order Altarelli-Parisi 



equations [25|. Since the BFKL evolution deals with an unintegrated gluon distribution, 



we calculate its derivative using the well-known singlet Altarelli-Parisi equation given by 

;*p(,. Q.) ^^['f [f) s(v. Q') + Pi? [f) nv. Q')) . 

where E denotes the quark singlet part, and Pgg\ Pgq^ the usual LO splitting functions. 
In the same way we produce a leading order MRS D_-type gluon, based on an input given 
in |6|. 



We also use a dynamically generated gluon distribution, for definiteness the GRV '92 
LO parametrization [^, that has the advantage of (a) being based on an explicit LO 
calculation and (b) being positive definite down to a low value of Q^. 



Furthermore, we take a look at the MRSA-Low gluon [^, which extends the valid 
range down to 0.625 GeV^ using an ad hoc form-factor-like ansatz similar to the one 
we employ for /. However, as MRSA is a NLO analysis, there is no consistent way to 
implement it in our BFKL evolution, which is neither an MS- nor DIS-renormalization 
scheme calculation. 



2.1 Treatment of the infrared region 

As mentioned in the Introduction, we need a suitable description of the infrared region; for 
our calculations we use the ansatz explained in detail in Jll], |l2l. This ansatz introduces 
three parameters: 
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• An IR- "cutoff" k^, i.e. a parameter which separates the infrared region, where an 
assumption on the fc^-behaviour of / has to be made, from the region where the 
Lipatov equation is solved numerically. 

• A parameter k"^ that controls the infrared behaviour of /. For k"^ < k1 we set: 

fi^^k^) = C'^^J{xX)- (4) 
The proportionality constant is given by 

k"^ + k'^ 

to guarantee continuity at k'^ = kl- This ansatz ensures that for /c^ — > 0, fi^x^k"^) ~ 
A;^, as required by gauge invariance. 

• A scale kl where we "freeze" the running coupling constant, i.e. 

as(fc^) — >as{k'^ + kl). 

This procedure applies also to our boundary condition /(a;o, k"^)-, although there is a 
slight modification to the pure AP gluon in order to soften its low k"^ behaviour: 

/^P(xo,fc^) — rP(xo,P + fcD- (5) 

The only purpose of the additional parameter kl is to ensure that we do not approach 
too closely a region where the gluon is unreliable; e.g. the DQ-typ^ gluon already ap- 
proaches zero at = 1 GeV^, and the Z}_-type gluon is not defined at all for such a low 
scale. 

Thus, 

This provides an infrared behaviour according to @, and for k"^ large enough /(xo,A;^) 
approaches /^^(xo,A;^). 
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In ref. [|TT|, the parameter is set equal to k^, whereas in |T2[ k'^ equals kf. We prefer 
setting k^ = k^ for reasons given below. 



2.2 Fixing kl 

It seems to be clear that the parameters introduced above should be small; so the simplest 



choice for these would be 12 



K = K = k'i = k'i = l GeV\ (7) 

However, there is a self-consistency constraint on the choice of k^, depending on k^ and 
Xq. The boundary condition @ should inversely be related to the (integrated) Altarelli- 
Parisi gluon g in Eq. ([Ta|) , that is 



k^ 



xog{xo,Q^)- I —f{xQ,k'^)=0. 



In a stricter approach, we apply the shift introduced in (^) also to the integrated gluon g 
in this equation. Then it is not very difficult to show that 



= xog{xo,Q'^+kl)-xog{xo,kl+kl). 

(9) 



It is obvious that this equation is fulfilled by setting k"^ = fc^, hence motivating our 
previous assumption. Using this and the asymptotic behaviour of g, 

xg{x,Q^ + kl) ^'xg{x,Q^), (10) 

we may return to (P) for simplicity, which gives us upon inserting our boundary condition 

(S) 

Jq /t -t- /t^ /t -t- /t^ 

k'^ + fc^ 

~ f^'' {xo^kl + kl)\Yi^^-^ + xog{xo,Q^) - x^g^xo.kl + kl). 



Set 


Agio [MeV] 




[GeV^] 




Do-type 


173.2 


0.01 


1.0 


0.95 


-D--type 


230.4 


0.01 


1.0 


1.51 


GRV '92-LO 


200.0 


0.01 


1.0 


0.19 


MRSA-Low 


230.0 (MS) 


0.01 


1.0 


0.44 



Table 1 : Parameters for various sets of parton distributions 

Thus we obtain an implicit equation for k'^: 

kl + kl ^ xog{xo,kl + kl) 

kl f^p{xo,ki + kiy ^ ' 

One can see that the value of fc^ mainly depends on the ratio g/f in the region of low 
Q^, and as the variation of / for different parton distributions is relatively small in that 
region compared to the variation of g, we conclude that it is basically the absolute value 
of g that determines the size of fc^. Hence, as a rule of thumb, the higher the value of 
xq g{xQ,Q'^) for small Q^, the smaller the resulting k^. It should be emphasized that, 
although the shift in Eq. (^) is taken into account in Eq. (j^) for both / and g, but in Eq. 
(pT|) only for /, the solution of Eq. (|n|) provides a good estimate on k^, very close to the 
value we get by minimizing the l.h.s. of Eq. for > 100 GeV^. The result for all the 
parton distributions under consideration is given in Table |l|, together with the value of 
A-QCD for four flavors used with each distribution. Figure |^ shows the unmodifled gluons 
and the modifled boundary conditions according to Table 1^. 

Regarding the Do-tjpe gluon, we see that the optimized value of k^ = 0.95 GeV^ 
is indeed very close to the naive estimate of 1 GeV^, as was already noticed in |12 . 



It was also mentioned there that a reasonable choice of k^ should lie in the range of 
0.5 — 2 GeV^, and we see that the D^-tjpe value lies well within this range, while the 
MRSA set is already close to its lower edge. The most extreme boundary condition in 
this respect is derived from the GRV parametrization, with a value of only 0.19 GeV^. 

As the more recent parton distributions favor a larger gluon g, the assumption of 
kl = I GeV^ (a good choice for the relatively small Do-type gluon) does not appear to be 
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the best choice for these gluons. 



Let us now emphasize the importance of the constraint (P) on k"^: Starting with the 
simple assumption (0) that all the parameters introduced should be equal to 1 GeV^, one 
gets a slope 

X= L,^I, . =0.5-0.6 



fd\n{l/x) 



x=10' 



for the unintegrated, BFKL evolved gluon distribution, as expected, no matter if Dq- 
type, D_-type, or GRV is chosen as input for the BFKL evolution. If we deviate each 
of the infrared parameters from (J^), we see that it is which has the biggest impact 
on A. With k^ decreasing, the slope rises, resulting in a slope A ~ 0.9 for GRV partons 
(fc^ = 0.19 GeV^) as the extreme limit. 

Keeping these considerations in mind, we will now concentrate our analysis on the 
structure functions. We will demonstrate that the slopes of / and F2 are related in such a 
way that for large A our calculated F2 is too steep to match it to the recent HERA data. 



3 Varying the parameters 



We construct the boundary condition at xq = 10 ^. The consistency constraint (^) 
determines fc^ and as a standard value for the other two IR parameters we choose kf = 
kl = I GeV^ The UV cutoff is set to 10"^ GeV^ and we stay with the Aqcd of the AP 
partons. We have written an evolution program to solve iteratively. Below k^ the 
necessary integration can easily be done analytically, above we use the GauB-Legendre 
quadrature. This transforms the integro-differential equation (|^) into a set of coupled 
differential equations, allowing us to use the standard Runge-Kutta method to calculate 
the evolution. The evolved gluon below xq is combined with the unintegrated AP gluon 
above Xq to obtain predictions for Ft and F^ by performing the integration in Eq. (|^) with 
Monte-Carlo methods. For convenience a fit of F2 = Ft + Fl is given in the Appendix. 
Finally we fit the background to the data, as discussed at the end of Section |0|, and 
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obtain our BFKL prediction. 



It is vital to check the dependence of the results on the parameters used for the IR and 
UV treatment. In Fig. |^ we varied all relevant parameters, using the Z^o-type gluon. All 
curves have been calculated at = 15 GeV^ and are already fitted to the shown HERA 
data with the soft pomeron background C/pX~°'°^. 

The strong dependence on k"^, which we expect from the corresponding variation in 
the gluon slope, is obvious. The curves for low are much too steep. At k^ = 0.2 GeV^ 
for example, the pure BFKL prediction is too high, even without any background. A 
background fit would then give a negative contribution, i. e. Cjp < 0, which is unphysical. 
In all such cases we set the background to zero. It is crucial, that k"^ can be precisely 
determined from the consistency constraint. Without the constraint, we could vary the 
slope of F2 from 0.5 to 0.8 by choosing fc^ within the shown range of 0.2 to 2.0 GeV^, 
rendering any serious prediction impossible. In the already mentioned comparison of a 
BFKL calculation with HERA data |]23[, k^ is treated as a free parameter. If the fitted 
kl should not be close to the consistent value by chance, it is doubtful that any strong 
conclusions can be drawn from a successful description of the data. 

The influence of the IR cutoff k'^ is comparably small. It should be kept in mind that 
kl has to be fitted separately for each fc^. Actually the induced variations in k^ dampen 
the dependence of the predictions on k"^ slightly. It is no surprise, that variations of the 
UV cutoff do not introduce much uncertainty into the predictions, as the running coupling 



already serves as an effective UV cutoff ||TT| , |T2|. The remaining free IR parameter, k^, 



has a sizeable effect on the curves. The slope of F2 varies from 0.53 to 0.58 in the shown 
range. The standard value of 1 GeV^ gives a slope of about 0.55, and the variations 
of kl represent an effective error-band of our calculations. 

A considerable dependence on Aqcd is expected and can be seen in Fig. ^. Fortunately 
this parameter is fixed, since we are using AP gluons with a given Aqcd as boundary 
condition. It is interesting to note, that higher values of Aqcd than the rather low 173.2 
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MeV used in the Z^o-type partons give wnfavourably steep slopes! Finally we take a look 
at variations of xq, noting that we have to determine separately again. We expect the 
steep curve for xq = 5 ■ 10~^ due to the long evolution length in x. But it is daring to use 
the BFKL equations at such high values of x and the data do not support such a choice. 
We can also see that the comparably flat curve created by a short evolution in x starting 
from Xq = 10^^ does fit the data well. But this success can obviously not be claimed by 
the BFKL evolution, since the data do not extend far below that xq. In order to test if 
the BFKL equations can describe the current data, we have to choose a larger xq ~ 10^^. 
The typical steep BFKL slope simply needs enough evolution length to develop. Small 
variations to lower xq from the usually used 10~^ do not affect the calculations strongly. 
The standard values fl^, |l^ used for and xq give slopes slightly below and above the 



average expected from the variations, respectively. Thus this choice of parameters will 
give a sensible prediction while still allowing comparisons with earlier calculations |Tl|, [12 . 



The results presented in this sections show that the dependence of the calculations on 
the parameters is under control. This statement would be impossible, if k^ did not obey 
the consistency constraint (|^). 

3.1 Using other input distributions 

Besides our detailed analysis of the Z^o-type gluon, we also studied the effect of feeding 
different parton distributions into our evolution, namely GRV'92 LO and MRSA- 



Low Q ^7^. The results for F2 can be seen in Fig. |. 



Obviously by using GRV and MRSA gluons as inputs for the BFKL evolution, a 
structure function is generated which is far too steep for the data. Even though we have 
got, in principle, some freedom in fitting an appropriate background, this is useless here, 
since already the pure BFKL part of F2 is too high for the HERA data, so that we must 
set the background to zero. 
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The steepness of F2 is closely related to the fact that the proper k'^ is much smaller 
than 1 GeV^, resulting in a larger slope A of the BFKL gluon /. As we explained above, 
this is mainly an effect of the size of the integrated gluon input g. This suggests the 
conclusion that the whole BFKL procedure has only a chance to work with the older 
(smaller) gluons. Modern parton distributions imply small values of the crucial infrared 
parameter fc^, and since this has a strong effect on the calculated structure functions, it 
is likely that one does not succeed in matching these to present experimental data. 



4 Comparison with data 

In Fig. 1^ we compare our calculations with very low data from the Fermilab E665 



experiment [|13|. Preliminary data of the HERA ep collider |]l5l at low are also shown 
in this figure. The published 1993 HERA data for low to medium are presented 
in Fig. ^ together with our BFKL predictions. All graphs show the BFKL calculations 
based on the Z^o-type and the D_-type partons, as well as the latest dynamical NLO 
renormalization group predictions (GRV), as presented in Also shown is the soft 
pomeron background, which is included in the Z^o-type curve. A general feature of all 
figures is that the difference between the Z)o-type and the D_-type BFKL predictions is 
very small after fitting the background. For this reason we do not discuss them separately. 
We have checked, that the R = Fl/Ft values used to extract F2 from the experimental 
data are close enough to those predicted by BFKL. Thus it is not necessary to reanalyze 
the data in terms of -Rbfkl- 



We first turn our attention to the E665 data [|T^]. The data do not extend very far 
into the small x region, making a check of the BFKL behaviour difficult. On the other 
hand the difference between the GRV (AP) and BFKL predictions is potentially large at 
small X. We also notice that the added background is comparable in size to the BFKL 
part at higher x. This flattens the steep BFKL behaviour, giving a good description of 
the data. But it is just the large contribution of the background which makes the very 
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procedure used doubtful. 

A further hint, that the successful description of the E665 data should not be taken 
too seriously, is provided by comparing the curves for the E665 data at 2.8 GeV^ and those 
for preliminary ZEUS data at 3.0 GeV^. While the E665 background is strong {Cjp = 
0.189) for Z^o-type, the optimal ZEUS background would be negative (C/p = —0.172) 
and therefore is set to zero. Thus the natural requirement, that the BFKL predictions 
for similar values should be approximately equal, is only fulfilled at very small x. 
For X > 10~^ the influence of the background quickly becomes stronger and the curves 
deviate. We conclude that BFKL can not describe both data sets consistently. It is also 
obvious from the figures that just this is possible using the usual NLO renormalization 
group equations, see the curves labeled GRV '94 in Fig. |^. 

Turning to the preliminary HERA data in Fig. ^ we find that the BFKL slope of F2 
is evidently too steep. The same tendency can also be found in the 1993 HERA data up 
to approximately 15 GeV^ as shown in Fig. ^ Everywhere in this region the background 
is very small or would even be negative, if it were allowed. Due to the larger spread in 
X, the preliminary data show the discrepancy rather clearly. The problem is rooted in 
the almost independent slope of BFKL, which is always > 0.5 and grows weakly with 
Q^. The dynamical GRV (AP) calculations give on the contrary fiat slopes at low Q^, 
which rise significantly with larger Q^. This is obviously in much better agreement with 
the data. For example, using the preliminary ZEUS data at 4.5 GeV^, we get a total 
of 15.3 (Do-type), 15.0 (D_-type) and of 2.81 (GRV '94 NLO) for the four data points. 

At higher than 15 GeV^, we see that the slope of the data becomes compatible 
with the one predicted by BFKL. The rising slope of the dynamical GRV prediction 
fits the data at least as well. It is interesting to note that even at higher a slight 
extension in the x range could provide an indication which evolution should be used. If 
future data should conform to the already visible tendency that a Lipatov-like slope is 
only obtainable at high Q^, say > 50 GeV^, then the simple LO BFKL formalism would 
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become implausible, since its validity at such high Q"^ is questionable. 

5 Conclusions 

We have shown that the AKMS method for calculating BFKL predictions of F2 remains 
stable under variations of the introduced parameters, if the consistency constraint Eq. 
(§) on fc^ is applied. As boundary conditions one has to choose older, i. e. smaller, gluon 
densities, since the large recent AP gluons lead to small /c^, which in turn produces overly 
steep slopes of F2. The BFKL boost at small x is simply too large for gluons constructed 
to produce the measured large F2 slope via the conventional renormalization group (AP) 
evolution equations. 

Thus we use the old Do-type and D_-type gluons with consistently fixed as input 
for the BFKL evolution. A further complication is introduced by the necessity to add 
a background contribution to the BFKL prediction for F2. Then it seems reasonable to 
require that the main growth of F2 is not driven by the chosen background, and that this 
background is comparably small. The last condition is not fulfilled in the region of the 
E665 experiment, casting serious doubts on the good agreement with the data. 

In the HERA region we find in contrast small background contributions, allowing for 
reliable comparisons with experiment. Especially the preliminary HERA data in Fig. ^ 
show that the almost independent BFKL slope is too steep for ^ 15 GeV^. But 
even up to the expected limit of applicability in of the BFKL evolution we find that 
the predicted slope is somewhat too steep. With improved statistics and maybe a slightly 
extended coverage in x, HERA should be able to assess the LO BFKL predictions for 
F2. If the tendency visible in the current data is an indication for future developments, 
we expect LO BFKL to fail the test. It remains to be seen whether extensions of LO 
BFKL — inclusion of the quark sector, NLO BFKL or theoretically consistent (energy- 
momentum conservation, etc. ) unified evolution equations — can improve the agreement 
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with the data. Our results also indicate that conventional (dynamical) renormalization 
group evolutions are still the best method for calculating and analyzing F2. 
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Appendix: Parametrization 



It should be convenient to have a simple parametrization of our theoretical results for 
F2. The following is a parametrization of the BFKL part only, to which an appropriate 
background still has to be added. 

We use an ansatz of the form: 

F^^^^ = ax-^ + l3x + -^, (A.l) 

which describes all curves shown in Figs. |] - ^ well. It is even possible to parametrize 
all pure BFKL results for the -Do-type and i5_-type gluons shown in Figs. ^ and ^ by 
choosing the following Q"^ dependence of the coefficients [t = \n{Q'^ / Aq^^))] 

a{Q^) = A + Bt + Ct^ + Dt\ 

f3{Q^) = Ef + Ft^, (A.2) 
7(Q') = G + H{QyGeY^) + It, 

A(g2) = j + Kt. 
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Do-type 


L>_-type 


A 


4.818 • 10-3 


8.07-10-^ 


B 


-2.460 ■ 10-3 


4.1 - 10-^ 


C 


7.386 ■ 10-^ 


1.579 - 10-^ 


D 


3.32 ■ 10-6 


4.370 - 10-^ 


E 


-2.23 ■ 10-1 


-1.89-10-1 


F 


1.15 • 10-2 


1.86 - 10-2 


G 


-5.3- 10-3 


-1.46 - 10-2 


H 


-6.66 ■ 10-^ 


-8.72 - 10-^ 


I 


6.54-10-3 


1.026 - 10-2 


J 


-5.113- 10-1 


-5.700- 10-1 


K 


-7.10- 10-3 


-2.89 - 10-3 


Aqcd 


173.2 MeV 


230.4 MeV 



Table 2: Coefficients defined in ([A.2|) of tlie parametrization (|A.1|) 



Tlie corresponding coefficients are given in Table ^ It is interesting to note that a small 
growth of A with has to be taken into account. 

It has to be stressed, that this parametrization is only valid within the range 10-^ < 
X < 10-2 and 0.8 GeV2 < Q"^ < 120 GeV2. On the one hand BFKL is not expected to be 
applicable even at the edges of this region. On the other hand we note, that the form of 
the ansatz has been tailored for this region only. The term (3 x, for example, will lead to 
wrong results for x > 10-2. 
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Figure Captions 



Fig. |l] Gluon ladder with tlie quark box for deep-inelastic ep scattering attached. The 
momenta are shown on the left side of the Feynman diagram and the parts of the 
/cT-factorization theorem, Eq. (|^), on the right. 

Fig. |2| AP inputs as indicated and the corresponding boundary conditions constructed 
via Eq. (^. The unmodified AP gluons are shown down to A;^ = 1 GeV^, whereas 
the IR treated inputs are continued below k"^. 

Fig. ^ Variations of the parameters. All curves are calculated at = 15 GeV^ with 
the Do-type gluon input. The background Cipx~^'^^ has already been added to the 
BFKL curves and statistical and systematic errors of the data have been added in 
quadrature in this and all following figures. The legend for is the same also for 
k'^ and k^. 



Fig. g Comparison of BFKL results, using MRSA-Low and GRV '92 LO gluon inputs, 
compared with data at several Q^. The corresponding curves for the Z^o-type gluon 
input is shown for comparison. 

Fig. |5| Low data of the Fermilab E665 experiment ||13| and preliminary data of the 



HERA experiments in comparison with our BFKL results. The background 
included in the Do-type curves is separately displayed. Note that GRV '94 refers 
to the conventional dynamical results based on NLO AP evolutions of valence-like 
input parton densities 



Fig. g As in Fig. | but using the HERA data of ref . . 
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